Unzipping Dynamics of Long DNAs 
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, The two strands of the DNA double hehx can be 'unzipped' by apphcation of « 15 pN force. We 

' analyze the dynamics of unzipping and rezipping, for the case where the molecule ends are separated 

I and re-approached at constant velocity. For unzipping of 50 kilobase DNAs at less than about 1000 

CNj ■ bases per second, thermal equilibrium-based theory applies. However, for higher unzipping velocities, 

,— ( ' rotational viscous drag creates a buildup of elastic torque to levels above fcsT in the dsDNA region, 

^ , causing the unzipping force to be well above or well below the equilibrium unzipping force during 

' ' ■ respectively unzipping and rezipping, in accord with recent experimental results of Thomen et al. 

' [Phys. Rev. Lett. 88, 248102 (2002)]. Our analysis includes the effect of sequence on unzipping 

CN ! and rezipping, and the transient delay in buildup of the unzipping force due to the approach to the 

steady state. 

■ 

55 ' I. INTRODUCTION 

^ . 

^ \ Double stranded DNA (dsDNA, the DNA double helix) is the genetic memory element of all cells. Two copies of 
G ■ the genetic information are encoded into the two complementary-sequence strands, which are base-paired together 
through most of the cell cycle. However, the two strands must be completely separated during DNA replication, and 
r-| , partially separated during DNA transcription. In cells, the separation of DNA strands occurs via forces applied by 
Q ' DNA-processing machinery. Force-driven dsDNA 'unzipping' is therefore of direct biological relevance. 
O I A few groups have carried out single-molecule studies of DNA unzipping by force. Essevaz-Roulet, Bockelmann 
and Heslot have studied the « 15 pN forces encountered during unzipping of 50 kb A DNA. Variations of the 

^-H ] unzipping force with sequence were observed, which are over the range of about 10 to 20 pN. Other groups have 
^ ■ carried out similar experiments on unzipping of DNA and on RNA helix-loop structures observing similar 

0^ ' unzipping forces. 

A number of theoretical works |^-[l^] have addressed the equilibrium statistical mechanics of dsDNA unzipping, with 
particular emphasis on the effects of sequence. Unzipping driven by DNA torque ('DNA unwinding') has drawn much 
less attention in spite of elegant experiments |lj] and corroborating theory jl^. As a result, theoretical consideration 
' of the combined effects of force and torque on unzipping have only recently been discussed [ p2[ . Since DNA unzipping 
involves rotation of the remaining double-helical DNA (Fig. 1), one expects that rotational drag torque should produce 
^ ' a contribution to the force needed to unzip DNA. Experimentally, no dependence of force on the rate of unzipping 
Ci ' up to about 1000 base pairs per second (bp/s) has been observed. However, in recent experiments of Heslot an 
^ appreciable increase (up to 40%) in unzipping force was observed, at unzipping rates in the range of 10 kbp/s fl^ . 

1 ■ This paper presents theoretical analysis of the velocity dependence of the unzipping force, for large molecules 
where the kinetics is dominated by viscous effects. To do this we introduce a dynamic model of unzipping which 
combines the polymer stretching dynamics of the extending single-stranded DNA (ssDNA) regions, the rotation of 

^ ] the remaining double-stranded (dsDNA) region due to the unwinding generated by the unzipping JTzl^p^, and the 
kinetics of the translation of the 'fork' separating the ssDNA and dsDNA reg ions. Although some work has been 
done on the Langevin dynamics of DNA unzipped by constant tension [^, p^|jl^ , a number of open questions remain. 
These include the unzipping force for constant end-to-end displacement velocity, the effect of sequence, and the role 
$— j of rotation of the dsDNA. 

Below, kinetic equations for unzipping show that beyond a certain unzipping rate, the predictions of the theory of 
equilibrium unzippin g ce ase to apply. For a 50 kb dsDNA this critical rate is about 1000 bp/s, similar to the threshold 
seen experimentally [ |l6| . In our theory there is a buildup of elastic torque in the dsDNA due to the drag torque 
associated with the dsDNA rotation. The ssDNA stretching degrees of freedom and the fork region itself remain near 
to equilibrium at experimentally accessible unzipping rates. We show that the relation between the unzipping force 
and the elastic torque buildup during opening is simply determined by the equilibrium coexistence between closed 
and opened and stretched base pairs, described by the binding-unbinding force-torque 'phase diagram'. 

We first describe the experimental situation in Section II, and then in Section HI we review the equilibrium theory 
of DNA unzipping, using a mean-field approach. We discuss the effects of force applied to the ssDNA ends, and 
torque applied 'upstream' to the dsDNA for homogeneous and heterogeneous sequences. In Sec. IV we discuss the 
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relaxation of ssDNA stretching and dsDNA twisting, and we then present a dynamic model for the propagation of the 
ssDNA-dsDNA 'fork' region for homogeneous and heterogeneous sequences. We first present a simple theory where 
we assume that the dsDNA twist is in a steady state. The problem of rezipping of a dsDNA is also considered, and it 
is shown that for rapid retraction, a left-handed viscous torque delays recombination of the ssDNAs. Finally in Sec. 
V we analyze the DNA twist dynamics in order to understand the delay in force buildup observed experimentally at 
high unzipping rates fl^ . 



II. EXPERIMENTAL SETUP 

The experiment of Thomcn, Bockelmann and Heslot unzips a A dsDNA of 48502 base pairs (bp) in 10 mM 
phosphate buffer, 150 mM NaCl, Ph 7, at room temperature. The A dsDNA is attached to two dsDNA linker arms 
each of 7000 bp; one of these is anchored to a movable glass slide, and the other is attached to a w-micron-diameter 
silica bead which is held in a laser trap (Fig. 1). The displacement of the glass slide at a controlled velocity in the range 
1 to 20 /im/s forces the molecule to open. After the molecule is opened (typically after « 25 /im of displacement) the 
stage motion is inverted, allowing the molecule to rezip into dsDNA (reannealing). The force transmitted to the ends 
of the linkers during the experiment is measured using the position of the bead in the laser trap, which is of stiffness 
kopt ~ 0.25 pN/nm. 

The stiffnesses of the two dsDNA linkers and the two ssDNAs depend on force, but are known from other experiments 
on dsDNA and ssDNA. At forces near 15 pN the total stiffnesses of the two dsDNAs (together 14000 base pairs) 
is kds ~ 0.1 pN/nm; the stiffness of the two ssDNAs depends on the number of unzipped base pairs n, and is 
kss « 70/n pN/nm. Below, we will neglect the « 0.1 /xm shift of the bead in the laser trap. 



III. STATIC UNZIPPING 



The basic unzipping force follows from the simple argument that unzipping will be thermodynamically favorable 
when the free energy required for unzipping one base, 1 < g < AksT, is equal to the mechanical work done, f£, where 
£ Ki 1 nm is the projected length of ssDNA liberated during unzipping of one base pair. The resulting force f ^ g/£ 
is on the order of 10 pN as it has been observed experimentally. Below we describe static unzipping in more detail, 
focusing on the force and number of opened base pairs at fixed displacement. The equilibrium theory of unzipping 
describes experiments for velocities small enough to allow the system to stay in thermal equilibrium. 



A. Homogeneous sequence 



We start by considering an idealized homogeneous sequence, with a uniform go — 2.5kBT, the averaged value on 
the A sequence. The free energy cost g per opened base pair may be obtained for the A sequence, for the experimental 
conditions described above, using the Mfold program [ [20| with stacking and pairing free energies measured by Santa 
Lucia |l|. 

To describe unzipping of the dsDNA we write down the work that must be done by the force to separate the ends 
of the linkers by a distance 2x. This includes the work done by the force to extend the two dsDNA linker arms by 
2xds, the free energy cost of opening n base pairs of the A DNA, and the work done by the force to keep the ends 
of the ssDNA regions separated by 2{x — Xds), and is therefore a function of the number n of opened base pairs, the 
extension 2xds of the dsDNA linker arms, and the total end-to-end distance 2x. We also include the work done by 
torque F applied to the end of the dsDNA region; this will be essential to considering the drag opposing rapid rotation 
of the dsDNA region. This free energy reads: 



Fx{n,Xds) = 2 Wds{xds) + 2Wssix - Xds,n) + n[go + TOq] 



(1) 



Here W{x) = f{x')dx' is the work done by the stretching force f{x) at fixed extension, for the dsDNA linker arms 
and the unzipped ssDNA, and 6*0 = 27r/10.5 — 0.60 is the number of radians of rotation made during opening of each 
base pair; this model is discussed in more detail in Ref. p^ . 

In the range of forces of 10-30 pN the dsDNA is extended enough that the leading contribution from its high-force 
entropic elasticity, plus linear stretching elasticity, give an accurate model: 



Xdsif) = Lds 



'2\fA ) ^ -Ids 



(2) 
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We use persistence length A = 48 nm and Young modulus jds = 1000 pN as determined by separate experiments p2[ . 
The length of each 7 kbp linker is Lds = 2.38 fj,m. 

The ssDNA is described by the freely-jointed-chain-like (FJCL) model ||2^, which gives the extension of n of the 
unzipped bases as: 

Xss = n Issif) (3) 

where 

fb \ keT 



coth 



ksTj fb 



(4) 



The monomer length d = 0.56 nm, the segment length b = 1.4 nm, and the stretching elastic constant -fss = 800 pN 
are fitted from experimental data Eqs. (j^ and (Q) can be used to obtain free energies via integration by parts: 

W{x)=xfix)- / x{f')df'. (5) 



Given the W^s and the total half-extension x, the minimum of the total free energy with respect to linker arm 
extension Xds, dF^in, Xds) /dxds = 0, equilibrates ssDNA and linker tensions. Then, minimization of the free energy 
with respect to the number of opened base pairs n, dFx{n,Xds) / dn — 0, determines the equilibrium unzipping force 
/„, via: 

2wss{fu)-ye^ = go (6) 

where Wss{f) = hs{f') df. 

Eq. (^) is a first-order transition coexistence condition, stating that work done by the force and torque in opening 
a base pair equals the base pairing free energy. Fig. 2 shows the solution of (|^) plotted in the torque-force plane. 
Note that overwinding torque in the dsDNA (F > 0) increases the unzipping force. The size of dsDNA torque needed 
to appreciably shift up fu is go/8o ~ 4:.2kBT; a left-handed (negative, corresponding to dsDNA unwinding) torque of 
this magnitude makes unzipping occur for zero force, close to the unwinding torque inferred from experiments [Q. 
Below wc will calculate how the overwinding transiently built up in the dsDNA during rapid unzipping, resulting 
from rotational friction upstream of the unzipping 'fork', will boost the ssDNA tension. 

The number of unzipped base pairs is simply calculated from the condition that the total displacement 2x is equal 
to the extension of the two linkers arms 2xds plus the extension of the two single unzipped single strands 2xss' 

X ^ nlssif) + Xdsif) (7) 

At the beginning of an unzipping experiment, the two linker arms first stretch until the extension Xdsu = Xdsifu) is 
reached. At this point the A DNA starts to unzip, and the force stays pegged at /„, with the number of opened base 
pairs proportional to further displacement: 

/ N ^ X- Xdsifu) , . 

^ Issifu) ^ ' 

The equilibrium force and number of unzipped base pairs for this homogeneous model, with go — 2.5kBT and at 
zero torque, are plotted in Fig. 3. The critical unzipping force is fu ~ 16 pN, and the average projection of each 
unzipped base along the unzipping direction is Issifu) = 0.44 nm. This is close to what is observed experimentally 
at small displacement velocities v < 1 /im/s; as unzipping proceeds the ssDNA is stretched to about 50% of its 
total contour length of « 1 nm/base. 

Given accurate knowledge of the elasticity of ssDNA, unzipping experiments can determine the pairing free energy 
at room temperature. For the homogeneous model, we find go — 2.5fcsT. Prior to these experiments, this free energy 
difference was indirectly inferred from model free energies obtained from study of DNA melting at temperatures of 20 
to 40 degrees above room temperature . A sequence-averaged point of view as presented above can give a rough 
account of thermodynamics of unzipping of large molecules. However, there is appreciable sequence-dependence of 
the base-pairing free energy. Rief et al have found that pure AT (the most loosely bound base pairs) sequences unzip 
at about 9 pN, while pure GC sequences (the most tightly bound base pairs) open at about 20 pN ||]. This range of 
force corresponds to base-pairing free energies of 0.8 to 3.8 fc^T per base pair. Therefore, a more detailed analysis 
of unzipping, and especially describing unzipping of short inhomogeneous sequences, requires models which take into 
account sequence-dependence and the cooperativity of strand separation p4] . 
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B. Heterogeneous sequence 



Sequence effects can be added by making go a function of n. The equilibrium opening of A-DNA at zero torque 
has been theoretically analyzed by Bockelmann, Essevaz-Roulet and Heslot ||^ , who have numerically calculated the 
thermal average of the force and the number of opened base pairs. They included thermal fluctuations of the ssDNA 
and dsDNA regions using a free energy of the form (P plus trap/cantilever elastic energy. Here we show how to 
obtain essentially the same results, using a preaveraging of go{n). 

The approach of the previous section computes the free energy of a given number n of opened base pairs, at 
fixed displacement x, using (|^) and (0). We estimate the fluctuations of the ssDNA, the dsDNA and the laser 
trap/cantilever, at « 15 pN using their combined stiffness ktot = [l/^ss + l/^ds + l/^opt]"^- This net stiffness 
decreases with the opening because of the inverse proportionality of kss to the number of opened base pairs. During 
opening of the first 5 kbp of A, the extension fluctuations in the length are « 20 nm, corresponding to 20 base pairs. 

To account for these fluctuations we have Gaussian-preaveraged the denaturation free energy g{n) using a standard 
deviation of 10 base pairs. To each configuration of n opened base pairs is associated a Boltzmann factor using free 
energy where the force is determined by the condition The thermal averages of the number of unzipped base 
pairs and of the force, as a function of the displacement 2x, are plotted in Fig. 4, for displacement after opening up 
to 5 /xm (« 5 kbp). 

The result of this calculation is in good agreement with experiment and the computation of Bockelmann et al. [p5| . 
The average slope of the opening curve gives the average extension of each unzipped segment along the unzipping 
direction 2Zss = 0.95 nm (one base pairs opens for each 0.95 nm of displacement). As in the experiment, the sequence 
generate a stick-slip motion; the opening fork stalls at G-C-rich parts of the sequence, giving a sawtooth pattern in 
the force signal and a step pattern in the number of opened base pairs [^SUlCj ]. 



IV. UNZIPPING DYNAMICS 



To describe the motion of the unzipping fork, we must consider four physical effects. First and second, we must 
consider the elongating dynamics of the dsDNA linker arms, and the unzipped ssDNA. Third, we must worry about 
the propagation of twist down the dsDNA; each 10.5 bases unzipped forces one more full right-handed twist into 
the unzipped dsDNA region. Finally, we must close the dynamical equations with a model for the translation of the 
ssDNA-dsDNA fork in n. For each of these processes we consider relaxational dynamics of the form: 

C- = -- (9) 
^ dt dz ^ ' 

Here z can be the position of a monomer of ssDNA or of dsDNA, the twist orientation 9 oi a, monomer of the dsDNA 
being opened, or the number of unzipped base pairs n. In each case F is the free energy, while ( is the relevant 
friction constant. For translational motion of ssDNA or dsDNA monomers, = Qn-qa; we take a — \ nm. For twist 
relaxation, Q = Airrir^ where r is the dsDNA hydrodynamic radius, and where = 0.001 Pa-s, the suitable value for 
aqueous buffer. We use r = 2 nm, twice the 'bare' chemical radius of the double helix; this generates force curves in 
accord with the experiment Note that r is the only adjustable parameter of our theory; it is essentially a friction 
constant. 



A. DNA stretch and twist relaxation times 



Since ssDNA and dsDNA are appreciably stretched by the > 10 pN forces applied during unzipping, their dynamics 
is reasonably described by local hydrodynamic friction. Logarithmic corrections due to long-range hydrodynamic 
coupling can be added to this discussion but without major effect. To estimate the order of magnitude of the dsDNA 
linkers and the ssDNA to equilibrate, we expand the stretching free energies around the force of / = 15 pN (the 
nonlinear elasticity of the previous section is used) to obtain: 

F^^Y. [^("' t)-z{n-l,t)~ z^f (10) 

n 

where z{n, t) is the position along the unzipping direction of either the nth monomer of the dsDNA, or of the unzipped 
ssDNA. In these two cases the monomer stiffness are either Kds = 1400 pN/nm (dsDNA) or Kss = 140 pN/nm 
(ssDNA) . Note that these should note be confused with the polymer stiffnesses kds and kss discussed previously. 
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The dsDNA twisting free energy is well described by the free energy of an elastic rod P6| in the range of torques 



relevant to unzipping experiments. This free energy may be written as in (10) where the degree of freedom z(n, is 
the twist angle of base pair n. The twist stiffness is = ksTC / /^"^ where the base rise A = 0.34 nm converts base 
index n to dsDNA contour length (note ksTC = 80 ± 20 nm is the usual elastic-rod twist rigidity j26|) 
The longest relaxation time of (g) with an elastic free energy ( [lO| ) is 

where either M = N^g — 14000 is the number of base pairs in the two linkers DNA, the number of unzipped base 
pairs M = 2ti, or the number of still zipped base pairs M = N — n. Substituting the relevant stiffness k and drag ^ 
(pi]), we obtain t^s = 3 • 10"'' s for dsDNA stretch relaxation, t^s = (2?^)^ x 1.4 • lO^^'^ s for ssDNA stretch relaxation, 
ttyj — [N — n)^ X 2 • lO"-'^ s for dsDNA twist relaxation times. 

Equilibrium will be reached for stretching or twisting if the relevant relaxation time is less than the unzipping 
time 1 nmxi„ « n/v. For the maximum velocities we are considering, 2v — 20 /im/s, tu — n x 5 ■ 10"^ s. Since 
tss/tu ~ n/10^, ssDNA in unzipping experiments where 2v < 20 ^m/s will be equilibrated until about 10^ bp are 
unzipped. Therefore dsDNA and ssDNA stretching is at equilibrium in experiments on A-DNA (n < 5 x 10^). 

By contrast, twist relaxation cannot reach equilibrium at the start of unzipping; the relaxation time at rt = 
is io = 5 • 10^'^ s. Thus, we will now describe the fork dynamics, treating the ssDNA and dsDNA stretching in 
equilibrium. We first develop a theory of steady-state twisting; in Sec. ^ we analyze the approach to this steady 
state. 



B. Fork motion 



The fork position n{t) will change as a result of the imbalance of ssDNA tension, opening cost, and dsDNA torque. 
Unzipping experiments occur at less than 20 kb/s, slow enough that each base is opened on average slower than 
single-base opening-closing times which are less than a microsecond. A quasistatic model of fork motion is plausible, 
with fork velocity in proportion to the free energy change associated with unzipping of one base. 

We write the relaxation equation for the fork motion in the continuum limit for the number of opened base 
pairs: 



dn 1 

Tn 



dt kuT 



{2w,s{f)-9o{n)^re„) (12) 



The single-base relaxation time r„ should be on the order of the diffusion time for the « nm-long bases r„ w 
QTirja? /kBT « 10~* s. The force f{n,x) is again determined from (^, and is implicitly a function of both the 
displacement and the number of open bases. 

Equation ( p^ ) indicates that if the ssDNA tension is large, the fork moves to larger n. The fork is static {dn/dt = 0) 
for the equilibrium state (^). Each unzipped base forces the upstream dsDNA to rotate through 6*0 — 0.63 rad. If this 
fork rotation is sufficiently rapid, the viscous rotational drag along the N — n dsDNA base pairs which remain to be 
unzipped, will generate elastic torque. 

In this subsection we assume that the DNA to be unzipped has reached a stationary state rotating at a uniform 
angular velocity u. Since each opened base forces a rotation of the dsDNA region by an angle 0o, we have u = O^dn/ dt. 
The viscous torsional drag for the dsDNA, treated as a cylinder of cross-sectional radius r and length A(A^ — n), is: 

T{n) = Am-ir'^l^{N ~n) LO (13) 

Combining (13) and ( p^ we obtain the equation of motion for the fork position: 

dn _ 1 2wss{I) - 9o{n) 



dt ksT Tn + {N — n)Tr 

where 



(14) 



Tr = = 2 • 10-^ s (15) 

The time Tr is comparable to the value expected for t„ (both are viscous times at the nanometer scale), but since 



Tr appears in (14) magnified by a factor k, N relative to r„, the rotational dynamics will be rate-limiting in most 
experimental situations. Also note that Eq. (|lj), and are independent of the value of the twist elastic constant C. 
Numerical integration of (14]) gives force and torque during unzipping. The initial condition is that unzipping begins 



when the force in the linkers reaches /„, i.e. n{t — Xdsu/v) = 0. 
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C. Rezipping dynamics 



In the experiment of p^ , following unzipping of « 25000 base pairs the molecule is allowed to rezip ('reanneal' 
in the nomenclature of biochemistry), by reversing the direction of the pulling velocity. If velocity is made negative, 
(1^) describes this process. During rezipping, the dsDNA rotates in a right-handed sense, generating a left-handed 
drag torque on the molecule. 

As shown in the phase diagram of Fig. 2 a negative torque promotes helix opening, decreasing the unzipping force. 
Therefore, the force during rezipping at high velocities is lower than the equilibrium unzipping force /„; the force 
drops progressively during rezipping because the rotational drag ( |l3|) increases with dsDNA length. Force during 
retraction is calculated via integration of ([1^) (with v — w), starting from initial condition n = 25000. 



D. Analytical estimate of unzipping force dependence on velocity 



We can estimate the increase of steady-state unzipping force with velocity using (y 
(7o = 2.5kBT. Changing variables from t to x = vt, (114) can be rewritten as, 



dn 
dx 



1 



(fcsT) V 



"^WssU) - go 

Tn + {N - n)Tr 



for homogeneous sequence 



(16) 



where f{x,n) is the equilibrium tension given by (7). In the limit of slow unzipping, the square bracket of ( [iq ) goes 
to zero, giving the equilibrium unzipping force (1) and the relationship ^ between extension x and number n„ of 
unzipped base pairs, n„(a;) = {x ~ Xdsu) /hs{fu)- 

We now suppose that the velocity is small (in a sense made more precise in the following), and that the unzipping 
force / and number n of unzipped base pairs can be expanded to first order in v, i.e. f{x) = fu + fii^) v + O(w^) and 
n{x) ~ nu{x) + ni(x) v + 0{v^). Plugging this into and ( p^ permits us to obtain two coupled equations involving 
/i and rii, with the results 



ni{x) 



, -I- [A^ - n„(x)]rr 

2[lss{fuW 



^ss ifu) 



(17) 



where h' denotes the derivative of h with respect to its argument. 
Therefore, at low velocity, the unzipping force reads 



fix) = /. 1 + 



v*{x) 



where 



v*{x) 



2 fu ^ss i fu) 



kBT[{N - nu{x))i 



(18) 



(19) 



The function v*{x) is plotted in the inset of Fig. 5. It is a rapidly increasing function of displacement x, bounded from 
below by v*{x'^^^) « 20 /im/s for A-DNA (N = 5 x 10^ bp). Therefore, one can expect to observe a large increase in 
the initial unzipping force for velocities larger than a few fim/s. 

Fig. 5 shows this force vs. velocity behavior which should be an upper bound to forces observed during the 
unzipping of a A-DNA (dashed line). The theory indicates a 10% increase in unzipping force as v is increased to 
about 2 /im/s, comparable to the initial rate of increase recently observed fl^ (velocities reported by Thomen et al. 
correspond to 2w, see Fig. 1). 

Eqs. (p^ ) and (^9|) also show that as the molecule unzips and N — n goes down, the torsional drag on the dsDNA 
is reduced, and the force needed to keep the fork moving goes down. For low velocity, this force drop will be a nearly 
linear function oi N — n^. This is the signature that the twist transport dominates the fork retardation. Finally, we 
note that (|l^) can be used to estimate how the force depends on velocity during retraction by simply inverting the 
sign of the velocity. 
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E. Rezipping at zero tension 



For sufRcicntly fast retraction the reanncaling of the double hehx, the rate of which will be limited by the rotational 
drag on the dsDNA region, will not be able to keep up with the retraction. As a result the force should essentially 
drop to zero. This effect will be especially pronounced during the later stages of rezipping, since the rotational drag 
experienced by the dsDNA region increases as rezipping proceeds. 

The velocity for 'free' rezipping under zero tension can be easily estimated, using the equation of motion for 
the dsDNA-ssDNA fork: 



dn 



1 



90 



knT 



Taking initial condition n{t = 0) — no, we obtain 

n{t) = N 



TV 



N + 



(TV 



no 



(20) 



go 



1/2 



(21) 



Recall that r„/Tr = 5 is not a large number. Therefore, for long {N > 10"^) and initially totally open {no = N) 
molecules, the rezipping follows N — n ^ (2 ffo t/iksT r^))^/^, where the exponent is determined by the linear 
dependence of the total rotational drag on the size of the rezipped domain. 

The total time needed to rezip for a 'free' fork is therefore T « (2A^ — no)no Tr ksT j (2 jtq), which for ri(0) — 25000 
is w 0.75 s. Therefore, for half-unzipped A-DNA, retraction velocities of 1v > 50 fim/s are at essentially zero tension. 
Equations (13) and (20) lead to the dsDNA torque for zero-tension rezipping: 



r = 



go 

Oo 



1 



rniN-n{t)) 



(22) 



with n{t) given by (|^). The maximum torque that appears at the fork during zero- force rezipping is the critical 
torque for the opening at zero force go /do, in accord with the phase diagram of Fig. 2. 



F. Results for homogeneous sequence 

We now present numerical results for integration of ( |l^ for the homogeneous case go — 2.5^^7, rj — 10^^ Pa s, 
r = 2 nm, A = 0.34 nm, and 9o — 27r/10. The ssDNA elasticity FJCL parameters are b — 1.4 nm, d = 0.56 nm and 
7ss = 800 pN. The dsDNA linker elasticity parameters are L^s ~ 2.38/im, A ^ 41 nm and = 1300 pN. The fork 
and torsional relaxation times are t„ — 10~^ s,Tr = 2 10~^ s. 

Force as a function of displacement after opening is shown in Fig. 6 for velocities 2v of 4, 8, 16 and 20 fim/s. During 
unzipping one observes an initial force upswing as the dsDNA linkers are first tensed, followed by a force peak and a 
gradual force reduction due to the decrease of torsional drag, as unzipping proceeds. For 2v — 20^m/s the peak force 
is fmax — 23 pN, corresponding to a 7 pN unzipping force increase relative to the equilibrium value fu =15.7 pN; 
this is in good agreement with the increase of « 10 pN observed in the experiment of fl^ . On the other hand, the 
initial force increase observed experimentally is smoother than the theory of Fig. 6. Below we will show how sequence 
effects and twist relaxation dynamics reduce the initial rate of force increase. 

Fig. 6 also shows force during retraction. For 2v = —4 fim/s, there is already a noticeable force hysteresis relative to 
the 2v = -t-4 /im/s extension curve. For 2v > —8 /im/s, the ssDNA force approaches zero at the end of the retraction 
cycle. Similar 'hysteresis loops' were observed by Thomen et al. [p^ . 

Fig. 7 shows the DNA torque at the unzipping fork, during these unzipping-rezipping cycles. The unzipping torque 
reaches a peak coincident with the force peaks of Fig. 6. At 2v ~ 20 fim/s the maximum torque during unzipping 
is Tmax = 2.8kB T, while the maximum unwinding torque during rezipping is Tmax = S.SkBT, slightly smaller in 
absolute value than the zero-tension limit of F = —go/do = —4.2kBT . 

The force as a function of torque during opening and closing run at the maximal velocity of 2v = 20 /im/s are 
included in the phase diagram of Fig. 2. These curves follow the ssDNA-dsDNA equilibrium transition, indicating that 
the opening and the closing of the single base pairs takes place essentially at equilibrium. The relation between torque 
and force at the moving fork is, under likely experimental conditions, determined by the base-pairing interactions in 
the same way as at equilibrium. Points A,B,C,D in Fig. 6 and Fig. 7 are mapped to the coexistence curve in Fig. 2; 
for example, the force peak corresponds to point B in the phase diagram / = 23 pN, F = 2.SkBT. Turning this 
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around, this indicates that the experimental force-displacement curve can be used to infer the torque-displacement 
curve, using the equilibrium coexistence line of Fig. 2. 

The peak forces observed during unzipping are plotted in Fig. 5 (solid line), and match the approximation of the 
previous section at low forces. Fig. 5 also shows the slightly lower force occurring at the point where the molecule is 
50% unzipped (dot-dashed line). This force is reduced simply because at the half- unzipped point, there is less dsDNA 
remaining to provide torsional drag than at the peak force point (see Fig. 6). 



G. Results for heterogeneous sequence 

Fig. 8 shows the result s of numerical integration of using the preaveraged A-DNA pairing free energy go{n) 
as discussed in Sec. IIIB, The force as a function of the displacement for the velocities 4,8,16,20 /Ltm/s is in good 
agreement with experimental data, and reflects the sequence; note, for example the progressive increase in G-C 
percentage in the first 2000 base pairs, the decrease from the base pair 20000 to 24000 and the steeper increase from 
base pair from 24000 to 24500, all these feature are well reproduced in the experimental and theoretical curves. An 
interesting effect is that the fluctuations in the force due to the sequence are attenuated, especially during rezipping, 
at higher velocity. 



V. EFFECT OF TWIST RELAXATION 



The initial increase of the force of Fig. 8 is still faster than that observed experimentally. We now examine the effect 
of the initial twist relaxation dynamics, focusing on its influence on the force signal at the beginning of unzipping. The 
combined set of equations for the fork and twist comprises a moving-boundary-condition problem that is difficult to 
solve even by computation. In this section, we construct an approximate solution for the combined twist and opening 
dynamics, valid when the number of unzipped base pairs remains small with respect to N. 



A. Memory kernel for fork motion 

We rewrite @ and (|l^) as 

80 _ 

dt ~ dm^ ' ^ ' 

where m is the continuous base pair index, and where Ttw = 47r?7r^A^/C ~ 1.410^^^ s. Prior to unzipping, the 
dsDNA is relaxed, 

6l(m,0) = (0<m<7V). (24) 

The boundary conditions expresses that the in — N extremity of the molecule is free (zero applied torque) , while the 
location and unwinding of the other extremity depend on the number n{t) of unzipped base pairs at time t, 

^iN,t)^0 , (25) 
9{n{t),t) ^ -9on{t) . (26) 



We now introduce the Laplace transform of the twist, 9{m, p) ~ dt e 9{m, t) . The solution of (|2^) with boundary 
condition (^5|) and initial condition ( [24| ) reads 

e{m,p) ^ J{p) cosh[q{N ^m)] , (27) 

where = prtw and J{p) has to be determined to fulfill the remaining boundary condition ([2^). Defining the Laplace 
transform n{p) of n{t), ( [26| ) reads 



-eon{p)^ dte-P^e{n{t),t) ~e{Q,p) , (28) 
as long as n{t) <C N from (p7|). Eliminating J{p) from (p7|) and (psh, we obtain 
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0{m,p) = -6*0 n{p) 



cosh [q{N — to)] 
cosh \qN] 



(29) 



From (p9[), the derivative of the twist just upstream of the fork, d9ld'm{n{t), t), determines the fork torque. Inserting 
this into the Laplace transform of the equation of motion (|l^ ) for the number of unzipped base pairs yields 



T„ n{p) = 



p + aq ta,nh{qN) Jq 



(30) 



with a = 9qC / l^/Tn- Inverse Laplace-transforming (|3Q), we obtain a self-consistent integral equation for the number 
of unzipped base pairs, 



n{t) = / dt' G[t - t') 
Jo 



2wss{f) - 5o(«(<)) 



where the memory kernel G is defined through 



G(r) -^r^ 



(31) 



(32) 



Tn J -oo ^TJ"* p + aq teL'Dh{qN) 

Let us stress that the force / in ( ^l|) depends on n and t through (0). The poles of (||) are located on the real 
negative semi-axis, at pi = —{ye/N)'^/Ttw', yi is the root (unique) of tany = —b y such that \ye — iirl < Tr/2 {£ > 0), 
and b = aNrt^. Calculation of the residues is straightforward, giving 



G(r) = e(T) 



where 0(t) = 1 if t > 0, otherwise. Finally, 



1 



^ 1 



Ar2 



b + yf/b 



(33) 



(34) 



gives the I* elastic relaxation time of the dsDNA region. The longest relaxation time for A-DNA is ti = 3.5 10~ s. 

If the number of unzipped base pairs n{t) is small with respect to iV, then (|l]) is valid. For A-DNA this condition 
happens to be true over the time range < t < ti where ti is the longest relaxation mode of the double helix [see 
(^)]. This means that the force will reach its peak at about ri, and then will not vary significantly for later times 

Ti. The dsDNA rotation dynamics reaches a stationary regime, and at later times can be considered to be a rigid 
cylinder rotating at angular velocity lj = 9odn/dt. For times beyond ti, the twist-relaxation dynamics reduce to just 
the differential equation for the number of unzipped base pairs (|lj) . 

To analyze the dynamics for t < ti we solve ( |3l|) iteratively. We start from the equilibrium zero speed solution 
no(i) fx t. At step i, the number of unzipped base pairs as a function of time, ni{t), is inserted in the r.h.s. of ( |3l| ) 
and ni+i{t) is collected on the l.h.s. The iteration is repeated until convergence is obtained, which takes about 20 
iterations. We also solved the ordinary differential equation (|lj) using numerical integration routines, and observed 
that the two curves match accurately for times i > ti as expected. 



B. Results 



Fig. 9a shows force versus displacement including the twist relaxation dynamics for homogeneous sequence {go = 
2.5kBT). The only difference with the curve obtained without the twist propagation (shown for comparison for the 
velocity of 20 /im/s, dashed curve) is the much smoother initial force increase. Theory indicates a force increase 
spread over the first 2± 2.5 fim of displacement; after this extension the result converges to the one obtained without 
twist propagation. While we do observe an initial 'delay' of the force increase, the range over which theory predicts 
this effect is shorter than that observed experimentally (about 5 /im). 

We have also studied the dynamics of the opening fork including both twist propagation and the A sequence. In 
presence of a complex free energy landscape for go, the iterative scheme exposed above does not converge easily to 
the solution of ( pT| ) . We have therefore consider smoother landscapes through a preaveraging of the sequence over SN 
bases, with SN ranging from 1000 down to 100. Such values permit to reach numerical convergence and are sufficient 
to detect sequence-induced effects on the micron scale. The resulting force signal is shown in Fig. 9b for maximal 
velocity 2v — 20/im/s. The percentage of GC bases increases during the initial opening, and this effect spreads the 
initial force increase over the initial 5/im of extension. This trend is in qualitative agreement with experimental 
findings, though the calculated force is lower than the experimental value by a few pN. 



9 



VI. CONCLUSION 



We have presented a theory of DNA unzipping dynamics, for kifobase or fonger dsDNAs. We have shown that 
torsional drag built by the rotation of the double helix around its axis is the dominant frictional contribution in 
the opening of A-DNA molecules. Easily observable nonequilibrium effects for A-DNA are expected for unzipping 
velocities in excess of 2 /xm/s. The results of our theory are in agreement with recent observations by Thomen et al 
of a roughly 40% increase in unzipping force for A-DNA unzipping at 2v = 20 iim/s. 

We have neglected a few physical effects in the discussion above. First, we have not explicitly included effects of 
transport of the dsDNA base-pairs to the fork. This is particularly relevant to the experimental setup of Refs. |]l|,pHl6| 
where the fork moves relative to one of the ssDNA anchor points. This may introduce an additional translational 
contribution to the dsDNA drag (for a rod model, again proportional to N — n). However, an estimate made by 
Thomen et al suggests that this force should be small relative to the unzipping forces. In fact, at the lower 
fork velocities ~ 2 /im/s where nonequilibrium effects are observable, separate experiments and theory shows that a 
A-DNA coil should only be shghtly stretched Q . 

A second factor that we have ignored is the possible effect of dsDNA intrinsic bends. Nelson has recently argued 
that such bends should induce an orders-of-magnitude effective enhancement in the rotational drag coefficient [ p9[ . 
In the present experiment, it appears that this effect is nearly absent. Possibly, the roughly twofold enhancement of 
r over its 'bare' chemical value, needed to generate the observed unzipping force enhancement, is due to permanent 
bends or other structural inhomogeneities along the rotating double helix. 

A third, and potentially interesting effect is that once appreciable torque is built up in the dsDNA, there is the 
possibility that the upstream dsDNA may begin to writhe pq |. For a dsDNA under zero force, writhing (supercoiling) 
occurs when |r| > ksT |^^. Tension in excess of ksT/A « 0.1 pN {A = 50 nm is the dsDNA bending persistence 
length) pushes the writhing threshold up to |r| > {4kBTAfy/^. Writhing (chiral coiling) of the dsDNA could increase 
the effective r, even without the formation of plcctoncmic supercoils. Formation of plectonemes is straightforward if 
the dsDNA region starts as a random coil, since there will be near-crossings every few kb (every 5 to 6 persistence 
lengths) to act as plectoneme 'anchors'. Once plectonemes form, we expect a large enhancement in the effective 
friction for dsDNA rotation, and a large increase in unzipping force. 

It would be very interesting to see results for an experiment carried out in the geometry of Fig. 1. This might be 
done using two translated laser traps, which would allow much less perturbation of the dsDNA coil during unzipping. 

A natural question is raised by the absence of a noise term in eqn (^2|), preventing the system from probing the 
whole free energy landscape at equilibrium (see also ||l^. Section VII). This approximation, which makes easier 
for the system to be blocked in a local minimum (stick regime), is expected to be valid at large velocities i.e. when 
the landscape changes very fast, and times scales are too small to allow for barrier crossing between stick and slip 
states. To further test the validity of eqn (jl^), we have calculated the number of open base pairs n as a function of 
displacement 2x at extremely low velocity e.g. 10 nm/s for the A sequence. Results are in very good agreement with 
equilibrium predictions of Fig. 4; the only difference is that unstick jumps (Inset of Fig. 4) take sometimes place < 5 
nm after their equilibrium counterparts. 

Finally, we note that the intrinsic fork motion time t„ might be larger than the « lO^^s assumed in this paper 
using dimensional considerations. The activation barriers to opening of successive bases p],p^ might make t„ larger; 
a recent estimate based on the analysis of a RNA opening experiment at constant force |5|] gave r„ ~ 2 10~^s [ pO| . 
However, there will be little consequence of a larger t„ for the phenomena discussed in this paper, which occur on a 
much longer time scale. 
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2 X (t) = extension of 2n ssDNA bases + dsDNA linkers 




FIG. 1. DNA unzipping experiment considered in this paper. Connections are made to dsDNA linkers attached to the ssDNA 
ends, and are used to pull adjacent 3' and 5' ssDNA ends apart. The ends of the dsDNA linkers move at equal and opposite 
velocities of magnitude v. As the dsDNA is converted to separated ssDNAs, the helical turns of the dsDNA must be expelled, 
forcing the remaining dsDNA to be rotated once for each 10.5 bases which are unzipped. 
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FIG. 2. Phase diagram of a homogeneous DNA molecule under applied force and torque. The base pairing free energy is set 

to go = 2.5fcsT. The double helix conformation is thermodynamically preferred when the applied force is smaller than some 
torque-dependent critical value, e.g. fu = 15.7 pN at zero torque. Denaturation may be driven by torque only. In the absence 
of applied force, DNA opens when an underwinding torque larger (in modulus) than 4 ksT is applied. Points A,B,C,D refer 
to the force vs. displacement curve of Fig. 6. 
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FIG. 3. Force (pN,full curve) and number of unzipped base pairs (kbp, dotted curve) at equilibrium for a homogeneous DNA 

molecule, as a function of the displacement x after opening initiation. Prior to unzipping (a; < 0), the force vs. extension curve 
reflects the elastic behavior of linkers. As opening proceeds {x > 0), force is constant at fu = 15.7 pN, for base pairing energy 
go = 2.5/csT. The number of unzipped base pairs increases linearly with x, with a slope ~ 1 bp/nm. 
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FIG. 4. Force (pN, top curve) and number of unzipped base pairs (kbp, middle curve) at equilibrium for the A-DNA molecule, 
as a function of the displacement 2x after opening initiation. Bottom curve is the base pairing free energy vs. index of base 
pair, from a Gaussian average over 20 bp. The number of unzipped base pairs increases linearly with x, with characteristic 
stick-slip steps (see inset). 
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FIG. 5. Theoretical unzipping force as a function of pulling velocity for homogeneous DNA. Star indicates equilibrium 
unzipping force. Dashed curve shows the 'steady-state' approximate formula (^sj), solid curve shows the peak force encountered 
during integration of (^). The dashed and solid curves have the same dependence at low velocity, tending to the static 
unzipping force fu ~ 15.7 pN at zero velocity. Inset shows characteristic velocity v* entering steady-state formula is plotted as 
a function of displacement 2x. For velocities above «*(0) = 20 fim/s, the steady-state formula starts to be appreciably above 
the peak force obtained by integration of (^^. The dot-dashed curve shows the force at the half-unzipped point when about 
25 kbp are unzipped. 
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FIG. 6. Force vs. single-strand extension x — vt obtained from integration of (0), for homogeneous DNA plus 7 kb linkers 
as in Fig. 1. Outgoing (pulling) curves for 2v — 4, 8, 16 and 20 ^m/s are shown (solid curves, bottom to top). For these rates, 
successively higher unzipping forces are obtained. In all cases an initial force increase associated with pulling the linkers taut, 
is followed by a force peak, and then a slow force reduction. The slow and linear reduction of force with extension is due to 
progressively less dsDNA being left to provide rotational drag to oppose fork motion. Force during retraction following the 
extensions are shown for 2v = 4, 8, 16 and 20 ^m/s (dashed curves). For 2 and 4 fim/s, relatively small hysteresis loops occurs. 
However, at retraction at > 8 yum/s the hysteresis is larger due to larger rotational drag, the force drops near zero at the end 
of rezipping. 
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FIG. 7. Torque in the dsDNA region immediately adjacent to the fork, vs single-strand extension x — vt, obtained from 
integration of (t4), for the cases shown in Fig. 6. Solid curves indicate extension, dashed retraction. Progressively higher 
torques are obtained at successively higher pulling velocities. There is an appreciable peak in the torque during the early stages 
of unzipping, followed by a gradual torque decay as the remaining dsDNA provides progressively less rotational drag. This 
torque buildup is due to the rapid fork motion forcing the dsDNA region just upstream of the fork to be under overtwisting 
strain. The torques obtained during retraction (dashed curves) show that rewinding of the molecule generates left-handed 
elastic torque in the molecule; for retractions of 8 and 16 fim/s this torque exceeds the 'free' limit of go/&o ~ —AknT discussed 
in the text. 
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FIG. 8. Force vs. single-strand extension x — vt for A-DNA from integration of ([l^), for 2ii = 4, 8, 16 and 20 /im/s (outgoing 
curves: full line, retracting curves: dashed lines). Note that fluctuations of the force due to sequence effects decrease with 
increasing velocities. 
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FIG. 9. Force vs. single-strand extension x — vt with twist propagation obtained from resolution of (^l|). (a) Result for 
homogeneous DNA {go — 2.5kBT) and with 7 kb linkers. Outgoing (pulling) curves for 2v — 2, 4, 8, 16 and 20 fim/s are 
shown (solid curves, bottom to top). For these rates, successively higher unzipping forces are obtained. The initial force 
increase associated with pulling the linkers is much smoother than without twist propagation (shown for 2v = 20 fim/s). (b) 
Corresponding results for inhomogeneous sequence for 2v = 20/^m/s, compared to the homogeneous case (dotted line). Curves 
correspond to preaveraging of the sequence over 1000 (dot-dashed line), 500 (dashed), 250 (long-dashed) and 100 (solid) base 
pairs. 
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